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SYNOPSIS 

Hciri Vasudeo Sahasratuddhe , Ph.D., 

metrical Engineerir^ Department, Indian Institute of Technology ,KnnDnr 
Hay 1968 

Models for Comiuter-aided Design 
(Sensitivity and Optimization Models) 

There is increasing interest of late in extending the 
possibilities of computer-aided design for the solution of a 
large .class of problems. A large amount of tho present ?.ork in this 
area has however been confined to network design problems’. An 
excellent cross-section of such work is found in the special issue 
of the Proceedings of IEEE on computer-aided design. 

Intimately associated ’v^ith this problem of computer-aided 
design are the relative merits of on-line computation versus batch 
processing'. In eithor of these, tho theoretical studies in 
devoloping proper sinulation models have -ngaged the attention 
of many research wodesrs. 

One important theoretical study of this nature is sensitivity 
analysis 'vhich constitutes an important link between analysis and 
design. In this thesis, sensitivity models have b 6 gn_j.evelopod 
which are algorithmic in nature and hence can easily be incorporatad 
into computer formulaticn and solution.: 

These GopaiTuter models for sensitivity analysis constitute 
an extension of the state-space formulation of linear and non- 
linear net v'^orks. while revierang this necessary material, a 
contribution is also made for a rigorous inclusion of non-linear 
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components which have non -mo no tonic characteristics. The equations 
for such componertts do not satisfy an essential condition to 
satisfy the sufficiency conditions for the exist ome of solutions 
to a system of non-lineox algebraic equations, This difficulty is 
circumvented by considering such components as dependent drivers. 

Thus, restrictions in the form of terminal equations are shifted 
to restrictions in topological considerations. These latter 
conditions axe stated. 

If design concepts are to be based on this analysis, further 
study of the models is necessary. Thus, for example, higher-order 
sensitivity coefficients are required for the formulation of the 
minimization of sensitivity of a network. Models for ■ determining 
higher-order sensitivity coefficients have also been d^^velopod. 

For establishing design ob5ectivc:.s based on sensitivity, we 
need measures of sensitivity. These measures are defined both in 
the frequency-domain as well as in the time-domain. The models tha.t 
are developed heru are equally applicable for either frequancy- or 
tim-G-domain analyses and therefore lend thomsslvos for design 
purposes based on sensitivity measures. 

In stability analyses, the question of sensitivity appears 
quite naturally. For example, Lyapunov defines local stability 
in terras of the boundedness of the sensitivity of the response to 
changes in initial conditions. These aspects have been discussed 
hero quite explicitly. 

Determination of nth order aettsitlvityi coefficients for a linear 
network depends on the res onse of the given network and a knowledgo 
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of th-.. sensitivity coeffieionts • of order (n-1 ) . However, the basic 
data that is required for determining higher order 'sensitivity 
coefficients is just the information about the n.ety?ork graph and 
its component terminal equations-, the rest of the information 
required for the final solution can be determined from these. For 
a non-linear network, however, additional informaticn a,bout the 
eomrenent terminal equations is necessary for every increase in the 
order of the sensiti-vlty coefficie rt . 

The next important problem considered in this thesis is the 
optimization problem. Again, the thesis is that it would be 
advantageous to develop network models for this study so as to 
render the formlation of the optimization problem conducive to 
digital computer use. A proper jetv/ork interpretation has been 
given to the co-state equations of a network, thus preserving tho 
structural features of the probl’;..m, which are usuolly not considered 
in control tfcory. 

Having systematized the foriiiula,tion of co-state equations 
in terms of network models, it is shovn that in certain studies of 
the optimizatioji problem, the existing network analysis programmes 
can be used for the solution of the resulting two-point boundary 
valiie problem. It is also shown that sufficient conditions for 
optimality, such as the ones discussed by Bobbins, do not admit 
a network interpretation. 

Another type of optimization study which has interested 
research workers is the parameter optimization problem. A component- 
oriented fornulation of thi s problem is presented. 
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The Component -oriented formulation has distinct 
advantages over other formulation procedxires whore the parameters 
do not correspond to the individual components. In the prooedure 
given hare, the need for algebraic manipulation is minimized, and 
error analysis becomes quite straight-forward. 

One of the interesting theoretical questions in parameter- 
optimization studies concerns the search procedures ydthin the 
parameter space ; it would be interesting to know that part of the 
parameter space in Thich the search procedures are effective. 

This question has been parti, ally answered in the thesis. 

Leeds and Ifgron have conjectured that if in a network 
designed for minimal sensitivity the number of components is 
increased it is always possible to reduce the sensitivity. This 
conjecture is shown to be incorrect by means of a counter-example. 

This thesis uses network language exclusively.' Hovrever, 
the ansilysis is equally ap'plioable to all discrete physical 
systtms. The range of applicability is a«fiieved by meE.ns of 
graph-theoretic formulation procedures which are quite valid 
for other systems also. Exclusive use of the network langu-age is 
made merely for the sake of simplicity in prosontation. 



1 . Introduction 


There is increasing interest of late in extending the possibilities 
of computer-aided design for the solution of a large class of problems. 

A large amount of the present -vTork in this area has however been confined 
to network design problems. An excellent cross-section of such work is 
found in (16), Other examples of interest in this field are (2,14j19)* 
Initiirately associated vrith this problem of computer-aided 
design are the relative merits of on-line computation versus batch 
processing. In either of these, the theoretical studies in developing 
proper simulation models have exgaged the attention of many research 
workers ( 7f9j '12) . 

One important theoretical study of this nature is sensitivity 
analysis -which constitutes an important link between analysis and design. 
Early -workers in this area (for example, (17)) have mostly developed 
models which are analytical in nature, but which are not necessarily 
amenable for transformation into digital simulation models, 

In this thesis, sensitivity models have been developed vrhich 
are algorithmic in nature and hence can easily be incorporated into 
computer formulation and solution.' The general nature of these analyses 
is quite, similar to -the underlying procedures of problem-ori aited 
programmes such as (7>13) so that the latter' s range of usefulness could 
be enhsnaed when these additional features are incorporated. 
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I'hen network analysis programmes are not endowed with the 
ability for sensitivity analyses, several simulations of the same network 
configuration with slight changes in component values are necessary before 
one could arrive at an understanding about the sensitivity of the netTiork 
response with respect to component values. Network programmes '4^5.011 
incorporate the sensitivity models obviate the necessity for several 
simlations of the same network. 

We have so far dealt with analytical aspects of sensitivity. 
However, if design concepts are to be based on this analysis, further 
study of the models is necessary. Thus, for example, higher-order 
sensitivity coefficients are required for the formulation of the' 
minimization of sensitivity of a network as discussed in chapter 4* 

Models for determining higher order sensitivity coefficients have also, 
been developed. These models also constitute a generalization of the 
procedure for determining first-order sensitivity coefficients. 

For establishing design objectives based on sensitivity, we 
need measures of sensitivity. These can be defined either in 'bhe 
frequency domain or in the time-domain. The models that are developed 
here are. -equ ally applicable for either frequency-or time-domain analyses 
and therefore lend themselves for design purposes based on sensitivity 
measures. Chapter 3 deals with these questions. 

In stability analyses, the question of sensitivity apoears 
quite naturally^ For example, Lyapunov defines stability in terms of 
the boundedness of the sensi-tivity of the response to changes in initial 
conditions. These aspects have been di ssussed quite explicitly in c’lapter 
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In chapter 4> 'fche optimization problem is discussed. &^a5.n, 
the thesis is that it would be advantageous to develop netvrork models 
for this study so as to render the formulation of the optimization 
problem conducive to digital computer use. In some restricted probleriis, 
it is possible to use the initial -value programmes directly for the 
solution of the optimization problem. In some other cases, it may be 
possible to solve the two-point boundary-value problem of optimization 
by repeatedly carrying out initial -value analyses, \7ith this motivation, 
a proper network interpretation has been given to the co-state equations 
of a network, thus preserving the structural features of the prob lorn, 
which arejusually not considered in control theory (11,16). This additional 
feature can also be incorporated into the existing network analysis 
programmes , 

Another type of optimizaticn study which has interested 
research ?roiicers is the parameter optimization problem (1 , 14 ) • A 
component-le^jel formulation of this problem is presented in chapter 4 • 
There is scope for computerization of this formulation. The component- - 
level formulation has distinct advantages over other formulation 
procedures using indirect parameters which do not correspond to the 
individual components. The illustrative example presented shows how a 
minimum sensitivity network can be obtained by this formulation. 

This thesis has used network language exclusively. However, 
since the analysis is throughout based on graph theoretic concepts, it 
is egu ally applicable to all discrete physical systems. Exclusive 
use of the network laiguage is done merely for the sake of simplicity 
■ in presentation. 
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2. Computer Formulation and Solution of Linear and Konlinalr lTetwo 2 i:s 

2.1 Introduction •' The development of the main hody of this thesis, 
viz. sensitivity analyses and optimization problems, presupposes a 
knowledge of time-domain analysis of linear and nonlinear networks, 
liirthermore , since in these inquiries structural aspects are preserved 
state models based on graph- theoretic aspects are found essential. 
Consequently,, these latter models are briefly reviewed in this 
chapter. 

2.2 Linear Graphs : This section is devoted to some definitions a-nd 
theorems of graph theory. The theorems are stated here without proof. 
Their proofs can be found in (15). 

Definition 2.2.1 Edge ; A directed line segment whose ends are 
distinct is called an edge. An edge may meet other edges at either 
end, but may not intersect, or be intersected by, other edges. 

Definition 2.2.2 Vertex ; A common end-point of any number of edges 
is called a vertex. The number of edges meeting at a vertex is 
called the degree of the vertex. Trivially, an isolated point 
(which is the end point of no, or zero, edge) can be considered a 
vertex of degree zero . 

Definition 2.2.5 Graphs and Subgraphs ' A collection of edges and 
vertices such that no two edges have a point in common which is not 
a vertex is called a graph. A subcoll ecticm of the edges and 
vertices of a graph is called a subgraph. 

Definition 2,2 .4 Path ; A sequaice of edges which satisfies the 
following conditions is called a ’path. 
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(i) Each edge in, the sequence except the first has a vertex in 
common with the edge which precedes it in the sequence, and 
each edge except the last has a vertex in common with the 
edge Tihich follows it in the sequence . 

(ii) Ko edge or vertex appears repeatedly in the sequence. 

Definition 2 .2 .5 Initial, final, ani terminal vertices of a path : 

That vertex of the first edge which is not common to the second 
edge is called its initial vertex. Similarly, that vertex of the 
last edge which is not common to the previous edge is called its 
final vertex. The initial and final vertices of a path are 
collectively called its terminal vertices. 

Definition 2.2.6 Circuit s If the initial vertex of a path coincides 
with its final vertex the path is called a circuit. 

DefinitA^n 2 .2.7 Connected Graph : A graph is said to be 
connected if there exists, a path between any two of its vertices. 

A graph G lAhich has n^ subgraphs such that 

(a) no two of them have a vertex in common, and 

(b) their union is the graph G, 

is said to be in n parts. Those subgraphs are said to be parts of 
G. A connected graph is clearly in one part. 

Definition 2 .2 .8 Cutset j -^or a connected graph G, a subgraph C 
which has the properties 5 

(a) removal of G from G will leave G in two parts, 

(^) no proper subset of the edges in C has the property (a), 


is called a cutset. 
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Definition 2.2.9 Tree ; A connected subgraph of a connected graph G 
containing all its vertices and no circuits is called a tree of G. 

In a connected graph G of v vertices any tree contains exactly 
(v-l) edges. The edges of a tree are called its branche s while tin 
remaining edges of G are called chords . The collection of all chords 
is called the complement of the tree, or its co-tree . 

Definition 2.2.10 Fundamental circuit (f -circuit) and fundamental 
cutset (f-cutset). 

In connection with a connected graph G and its tree T, a 
circuit containing exactly one chord is called a f-circuit. Each 
chord belongs to a unique f-circuit. The number of f-circuits is 
therefore (e - v + l) where e is the number of edges in G. 

A cutset containing exactly one branch is called an f-cutset. 
Again, each branch belongs to a unique f-cutset, and hence the 
number of f-cutsets is (v - t) . 

Definition 2.2.11 : f-cutset and f-circuit matrices : Given a 
connected grajh G and a tree T, we now define two matrices A and 
D, called the f-cutset matrix and the f-circuit matrix respectively. 

The rows of A correspond to f-cutsets defined by the branches 
and its columns correspond to all edges. The typical element a. . is 
defined as follows s 

a^^ = +1 if the jth edge is in the f-cutset of the ith brancli^ 
and their ori oat ations coincide, 

= -1 if the jth edge is in the f-cutset of the ith branch 
and their orientations are opposite, 
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if the jth edge is note in the f-cutset of the ith 
"branch. 

Vvhen the columns of A corresponding to the tree are placed 
first and when the oid.er of these columns coincides with the row 
order, A has the form ; 


A = 


in 



2 .2 .1 


where U is a unit matrix of the apxropriate size*. 

The rows of B correspond to f-clrcuits defined by chords and 


its columns correspond to all edges, 
defined as follows : 


The typic al el ement b^ ^ is 


b. . = +1 if the jth edge is in the f-circuit of the ith 

chord and their orientations coincide, • 

= -1 if the jth edge is in the f-circuit of the 

ith chord, and their orientations are opposite, 
= 0 if the jth edge is noit in the f-circuit of the 

ith c hord. 

?Jhen the columns of B corresponding to the co-tree are 
placed last and vvhen the order of these columns coincides with 
the rocfi order, B has the fom : 



2 , 2.2 


Ijffcen the columns of the matrices A and B are matrices satisfy the 
orthogonality condition 

*The symbol IT will always be used in this sense. 
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ba' =0 2.2.5 

v/here the prime denotes transpose ; and, consequently, if A and B 
are in the form of Eqs. 2 . 2.1 and 2 . 2.2 we have 

= -^2 2.2.4 

2.3 Graphs of Networks ; ' 

A directed gdaph can be associated vdth an electrical netoXEck 
with lumped parameters. In this section \fe shall discuss this 
correlation. A detailed discussion can be found in ( 5 ) » 

Each edge in a network graph corresponds to a pair of 
voltage and current measureraaits between two terminals of some 
component (see Pigires 2 . 3.1 and 2.3.2). A two-terminal component 
is represented by a single edge, and in general an n-terminal 
component is represented by a tree (also called a "terminal graph") 
in tbje n-terminal complete graph. 

The cutset and circuit postulates governing the currents 
and voltages of the network are ; 

A I(t) =0 ■ 2.3.1 

and B V(t) =0 2.3.2 

v/here the matrices A and B refer to the network gra^, and I and Y 

are vectors ’jhose elements are the current and voltage variables 

associated liOith the edges of the netvTOrk graph*, arranged in the 

*Por simplicity of expression, these current and voltage variables 
will hencefourth be called, the currents and voltages of the graph. 
The distinction betv/een the mathematical model and its physical 
counterpart will not be blurred by such references made in the 
proper context. 




Birecti on-convention for 
corresponding edge in 
the network grajh . 


Birection-convention for 
volt .ogc and current 
measurement 

Pig. 2.5^1 Correlation between component 
network graph. 



Sufficient measaremt.nts for 
measuring transistor 
char actari st ics 


measurement and 



(b) 

Graph of a 
transistor 


Pig. 2.3.2 Graph of a maltiterminal component 
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same order as the columns of A and B. 

Equations 2.5*1 and 2.3*2 together provide e independent 
relations in 2e variahles; an equal number of independent relations 
in the same 2e variables is contributed by the component character- 
istics discussed in the following sections. 

2.4 Component Characteristics : 

First o^e shall consider a restricted class of components, 
namely, linear components. By describing a class of components v/e 
would also have described a class of networks - networks containing 
only comronents of the described class. Next, we shall move on to 
a more general set of nonlinear characteristics. These are defined 
in such a manner -that the linear characteristics become a special 
case of the corresponding nonlinear ones. A network of nonlinear 
components may therefore in general contain linear components also. 

The solution of lin-ar networks in both time-and frequency- 
domains is already textbook material (see, for example, (5))* A 
method of analysis for networks of nonlinear two-terminal components 
was given in ( 13 )* 

Furthermore, the formulation of networks containing multi- 
terminal aLgebraic components is also included here. The characteris- 
tics of such components are obtained as. an extension of the 
chaffiict eristics of two-terminal resistors. However, it is someti.n3s 
more convenient to replace raultiterminal algebraic components by 
the-ir equivalent networks o f two terminal components for the 
purposes of analysis. In fact, such representations are made use of 
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in compiter studies, shall make use of hoth these representations 
in our further study , 

2.5 Linear Components : 

Four types of tvo-terminal components and one type of 
fflultiterminal component come under this classification. They are; 


(a) Two-terminal resistors. 

The terminal equations can be written either in the 


form 

v(t) = R i(t) 2.5.1 

or in the form 


i(t)=Gv(t) 2.5.2 

(b) Inductors. 

The general form of their terminal equations is 


Y(t) 



2.5.3 


(c) Capacitors. 

The general form of their terminal equations is 


2.5.4 

In the above equations, the parameters R,G, L and C may 


bo positive constants or positive, bounded functions of time. A 
component whose parameter is constant is said to be time-invariant . 
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(d) Independent drivers. 

The voltage or the cuirrent of such a component is a 
specified function of time. The general forms of thifiir terminal 
equations are ^ 



i = i(t) 

2.5.5 

or 

V = ■v(t) 

2.5.6 


A driver whose voltage is specified is called a voltage 
driver, end a driver v/hosc current is specified is called a current 
driver. By definition, independent drivers are considered to he 
time-invariant components. 

(e) n-terminal resistors. 

• Such a component can he represented hy a tree grajhof 
(n - l) edges and the terminal equations 

X(t) = D Y(t) - 2.5.7 

corresponding to these edges. 

The notations X(t) and Y(t) for the column matrices 
(instead of 'V’(t) and l(t)) are used to allow for the possihility for 
the voltages and currents appearing in a mixed form. 

The vector X(t) in Eq. 2.5.7 consists of voltages of some 
(possibly none or all) of the edges and the currents of tie remaining 
edges 5 and Y(t) is the vector of the currents of those edges whose 
vol'.agfcs are in X(t), and the volteges of those whose currents aro 
in X(t). Such vectors X(t) and Y(t) are said to he complement ary . 
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In Eq. 2.5.7r^ is a matrix whose elements may he either 
constants or functions of time. If all the elements of B are 
constants, then the component is said to he time-invariant. 

A network of only linear components is said to he linear. 

In addition, if all its components are time -invariant then the 
netvork is said to he linear time-invariant. A network that is not 
time-invariant is said* to he time-varying. 

In another representation of a multiterminal component we 
replace it for the purposes of analysis, hy a network of two- 
terminal resistors and dependent drivers. Only two types of depen- 
dent drivers are necessary to he ahle to represent any raultiterminal 
component of the type described hy Eq. 2. 5.7* In- addition to those 
t’.vo more types of dependent drivers are defined for convenience. These 
four types are given hclow with their characteristics. 


Type 1 

V2(t) = ^ v^(t) 2.5.8 

Type 2 

V2(t) = ^ i^(t) 2.5.9 

Type 5 . ' ■ 

= H ' ■^1(1^) 2.5.10 

Type 4 

i 2 (t) - S i^(t) 2.5.11 


Onlj types 2 and 5 essential, since types 1 and 4 
be derived from them. 
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Example 2 ,5-1 

1'he raaltiterminal component repro scanted by the graph of 
Figare 2.5.1(a) and thJ terminal equations 

(t) 2.5.12 

can be replaced by the network of Figure 2.5.1(b) using dependent' 
drivers of types 1 and 4* Alternatively; it can also be represented 
by the network of Figure 2.5.1(c) when only dependent drivers of 
types 2 and 3 axe used. 

2.6 Nonlinear Networks 

The generalized terminal equations for nonlinear component 
chereoteristics are given below. Independent drivers, which are not 
mentioned belo’v, remain unchanged. 

(a) Two-terminal resistors 

2 .6 .l(a) 
2.6.1(b) 

2 . 6.2 


2.6.3 


or 


v(t) = R(i(t) ;t) 

i(t) = G(v(t);t) 
(b) Inductors 


<t) = 


_ M 


dt 


i(t) = T(jZi(t); t) 


(o) Capacitors 
i(t) 


M. 

dt 


Y 



1 

(t).= 

: 11 12 

^2 ! 
L- 


/ 2 I ^2_ 


v(t) = S(q(t) ;t) 




(^)) 

Network using dependent drivers of types 1 and 4. 
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(e) Multi terminal resistors 

X(t) = D(Y(t); t) 2 . 6.4 

Here, again, dependent drivers may be introduced to replace 
multiterminal resistors. The four-way classification of dependent 
drivers used in the linear case v;ill be followed here. These non- 
linear characteristics are represented by the equations ; 


Type 1 

■V2(t) = V(v^(t) ; t) 2.6.5 

Type 2 

■V2('t) = ■V(i^(t) 5t) ^ 2.6.6 

Type 5 

i2(t) = l(v^(t); t) 2.6.7 

Type 4 

i^Ct) = l(i^(t )5 t) ' 2 . 6.8 


When the function (R, T,S,D,Y, or l) in the charaoteri sties of a 
component does not depend on time, the component is said to be timo- 
invariant. ’//hen a network contains only time-invariant compone-nts, 
the network is said to bo time-invariant’. 

2.7 Overspecification, Degeneracy, and Restrictions of Structure. 

The existence of a circuit of voltage drivers or a cutset 
of current drivers amounts to overspecification. The existence of a 
circuit of voltage drivers and capacitors or a cutset of current 
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drivers and inductors is 'termed as a dSi^eneracy . A degeneracy may 
give rise to infinite currents or voltages in the solution of the 
nfi-twork. The following analysis presumed the absence of over- 
specifications and degeneracies. For such networks, tho -network 
graph will contain a tree "hich includes all the voltage drivers 
and capacitors and excludes all the current drivers and inductors. 

Such a tree is said to he mammal. 

The A and B matrices defined for a maximal tree are used ' 
for the formulation. This tree is called the formulation tree. The 
voltages of the branches and the currents of the chords of the 
formulation tree are called the primary variables while- the remain- 
ing voltages and currents are called the secondary variables. The 
secondary variables can be expressed in terms of the primary variables 
using Eqs. 2.5.1 and 2.3.2 as follows s 

^2^^^ = V^(t) 2.7,1 

and 

I^(t) = -A^2 ^2^^^ 2.7.2 

where tho subscripts 1 and 2 of V and I denote the tree- and co-tree - 
variables respectively. 

For the proofs of the results stated in this section sea 

(5). 

2.8 network equations for linear networks. 

when tho restrictions of tho last section are imposed on 
the network graph, we have six classifications of edges : 
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1. Branches corresponding to voltage drivers (dependent as ^*rell 
as independent drivers). 

2. Branches corresponding to cappcitors. 

5« Branches corresponding to resistors. 

4* Chords corresponding to resistors. 

5. Chords corresponding to inductors. 

6. Chords corresponding to current drivers (dependent as well 
as indopaidont drivers) . 

On the Basis of the six-v/ay classification, the V ani I 
vectors ind the pro multi plying A and B matrices are partitioned 
similar ljt*£ Correspondingly, VQ ro'wxLte the f-cutset and f-circuit 



and 



2 . 8.2 
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The terminal equations of capacitors and inductors can be 
grouped as follows : 


dt 




C'^ 0 


0 


-1 





2 



J 

5 


(t) 


2.8.5 


Usirg Eqs. 2.,8..1 and 2.8.2 we can rewrite Eq. 2.8.3 in terms 
of primary variables as follows 

r..-1 


d 

dtj 


C 0 

0 


0 -A, 


22 


"®22 ° 


V, 


(t) + 


0 -A, 


■®23 ° 


L^J 


(t) 


'"23 


-®21 ° 


(t) 


2.8,4 




Among the variables on the right of Sq. 2.8.4, ''?^2 ^5 

together constitute the state vector of th« network, namely, the 
capacitor voltages and inductor currents. 

The solution of Eq. 2.* .4 depends on the determination of 

. We shall discuss the 





rv,! 

the Yariahles in the vectors 

and 

5 


- ^ 6 ^ 



determination of these variables next 

•V 

f 

The elements of the vector 


1 


L ^6 J 


are the primary variables 

of drivers. If all the drivers in the network are independent, 
these variabl can be determined as functions of time. TiThen 
dependent drivers are present, however, the values of the variables 
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in the vector 




can he found if : 


(a) the dependent drivers depeni only on the primary 
variables of ot 1 b r drivers, capacitors, and indue tors j Cvnd 

(b) it is possible to arrange the dependent drivers in a 
sequence such that no dependent driver depends on any driver which 
follows it in the sequence. 

If these conditions are satisfied, the primary variables of 
the depaidait drivers can be determined in the order of any sequence 
viiich satisfies condition (b) above. 


trample 2,8.1 

Figures 2, 8.1 (a), 2.8.1(b) show networks in vdiich the 
conditions (a) and respectively, are violated. Figure 2.8.1(c) 

shows a netv©rk ih which these conditions are satisfied. Condition 
(13) is satisfied by the sequence (l,2). Once the voltages of tlB 
capacitors are knov/n, the voltage of driver 1 and the current of 
driver 2 can he found in that order; since we can irrlte 


v^ = o( v^ , 

and 

15 = ir(v5 - v^) 


v^ , v^ and v^ are all knom at the time when 


Next we shall consider tv/o procedures 

W, 

of the elements of the vector 


they are needed, 
for the determination 
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< 









(a) 



Fig. 2.8.1 Wet’^orks for demonstrating the determination of 
the primary variables of dependent drivers. 
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The resistor equations can be handled by vnritirg' their characteristics 
ezTlioitly in the primary variables : 


r 


r* 

L 



R, 0 
5 



-I 3 - 

1 

i 

1 



2 . 0.5 


Substituting the values of and in terms of primary variables 

f 

and transposing the term containing 


U + 


R, 0 


0, 


^4j 


to the left, wa get 


A. 


B 


15 


51 

0 


T, 


L^4 


~R, 

% 

f 

^0 -k ' 




0 -A.,., 



5 



52 


2 


55 


1 



1 




f + 




% 

G 

1 

0 

( 


Ir- 


-B. . 0 


I, 


4 J 

L L 

_ 12 J 


- 5j 

1 

- I"! J 


6 


2 .|..6 


It has been proved in (5) that the coefficient matrix on 

the left of Eq. 2.8.6 is nonsingular. Promult iplying Eq . 2.8.6 by 
the inverse of that matrix and substituting the resulting value of. 

into Eq. 2.8.4 'we get the dynamical equations (state equations) 



for the network ; 
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- 2.8.7 


Hhe re 



VvTien, dealing with nonlinear networks, the above formulation 
of resistor variables (Eq.s. 2 .8 .-5 to 2,8.8) is not possible, A 
formulation similar to tte nonlinear case is given below : 

First, the characteristics are v/ritten explicitly in the 
currmt s as folio v;s : 



2.8.9 


A part of the f -cut set equations (Eq. 2,5,1) can be re- 


written making use of Eq, 2.8,9 as ; 



2 . 8.10 
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Substituting the values of the secondary voltages 1 from 

4 

the f-circuit equations (Eq, 2.5.2) we get 


TJ 




0 ~ 

f 

0 ^ 


0 


” U 

♦ 

‘’J 

il 



-®12 





+ ^5 + “^53 ^6 = ~ 


2 . 8.11 


Once the solution to Eq. 2.8.11 is determined, we can 
solve the Network by integrating Eq. 2.8.4. 

We can treat Eq, 2,8.11 as a special case of the nonlinear 
resistor equations (Eq, 2 , 9 «. 5 ) of the following section, and use tho 
algorithm described there obtain its solutions. 

2.9 Extension o f Network Equations to Nonlinear Networks. 

The general nature of the formulation is the same as in the 
linear case. The equations of capacitors and inductors, having 
the forms of Eqs. 2.6.2 and 2.6.5 osn be grouped tegethea: as 
follows s 



'«2‘ 



d 




dt 





Substituting the value of ihe right-hand member of Eq. 2.9-1 
in terms of primary variables we get ; 
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L. 

dt 




0 -A 22 





' SS 




4- 



r®22 ° - 


- V 

1 ‘ 

w 


• .-A, 


'21 


"25 


L ^4- 


c~ ^ 

0 -A„_ 


Y. 

23 


1 

_ -B 2 , 0^ 


. " 6 . 


2 .9.2 


The elements of | 1 can. be found by the nonlinear part 

of the capacitor and inductor characteristics, having the forms of 
Eqs. 2.6.2 and 2.6.3, which can be grouped together as follows : 




r 

ro 

1 





1 

^ 5 - 


2 . 9.3 


It may be noted that the ith component of the function D, 


depends only on the ith element of 


0. 


5i 


The conditions on dependent drivers remain the same as 
in the previous section, and their characteristics can be simply 

used for obtaining the elements of 


1 


L^6 J 


. The new form of the 


tormina! equations of resistors is 




f 

f 

3 

= D 


3 


r 



I. 


\ 


u 4 J 


\ 

^4-1 


2.94 
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For the pose of this section, Eq., 2 . 9*4 contains only two -terminal 
resistor characteristics, and therefore the ith component of the 


function depends only on the ith comvonent of 



However, it 


is possible to group multiterminal resistors also into Eq, 2.9.4, in 
which case the above condition will not hold. 


From Eq. 2 . 9*4 we can obtain equations similar 


These equations are : 



/ 



- 0 - 


fij -1 

I u 

D I 


’'^0 + 


Y, -f 





2 

^ ! 

--^1- 

1 

— 1 


to 


5 


2.8,t1. 


■*"■^^2 ^5 ^55 “ .2 2.9*5 

Methods for solving Eq. 2 . 9*5 have been givfin by (5*4)* 
However, the conditions imposed by these authors are restrictive* 'They 
are s' 

For the typical resistor with a characteristics of the form 
of Eq. 2.6.1(b), the following should be satisfied ; 

(i) V2 > &(v2) 2.9.6 

(ii) I I ^ I '^1 “■ ^ 2 / 

for some positive, finite and M2 snd all v^ and 

These conditions can he stated in a short form hy saying tha,t 
each resistor should have a strictly monotonically increasing 
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characteristic with bounded slope. Since in such a case it is 
always possible to write the characteristics in the form of Eq. 

2.6*1 (b)* it is not necessary to consider the other form. The 
prorf for the above result can be found in (0). 

fharact eristics involving negative resistance, such as a 
tunnel diode characteristic, can not be handled by the algorithms 

of (3»4)- • 

Im such eases,, the characteristics can he t)!)ritten in only 
one of the two forms represented hy Eqs. 2.6.1(a), 2.6.1(h). For 
example,, the tunnel diode eharact eristics can only he written 
explicitly in the current variahle. Components which do not 
satisfy the conditions of Eqs,. 2. 9 *.6 , 2. .9 *7 handled if : 

(i) The variahle in which the characteristics is explicit 
can he mede the primary variahle^ and 

(ii) The secondary variable of the component can he 
determined from the knowledge of the primary variables of drivers, 
capacitors and inductors only. 

If the above conditions are satisfied,, the component can 
ho represented by a dependent driver controlled by its own secondary 
variahle*. Figure 2.9*^ shov/s an example of a network in which a 
tunnel diode is represented as a dependent current driver . 

2.10 Small-Signal Characteristics of Nonlinear Components : 

In chapter 4 we shall make use of the small-signal 
behaviour of nonlinear components. The small-signal behaviour of 
components is described in terms of the matrices D, and I 

, ■ ' .'Cl . T 



j 

— onfTTTr- — 

I 

+ JL 

-n 

r J 


(a) 


A network containing a tunnel diode 



l*ke network of witk tlie tunnel diode represented 
a dependent current driver. 

Fig. 2.9.1 A network containing a tunnel diode. 
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defined as follows ; 




2 . 10.1 




■ \ 

T^dioro X = 2.10.2 

As pointed out in the previous section, the ith component of 
depends only on the ith •omponent of X. Thtarefore, is diagonal . 

s 

The ma.trix is defined in terms of a more general form 
of the resistor characteri sties than Eq. 2.9»4- This form is 


X(t) - (Y(t)) 2.10.3 

where X and Y are complementary vectors of voltagvrs and currents, 
as explained earlier. From th^a^ the small-signal characteristics 

S' 

E is defined as : 
r 





2.10.4 


If Eq, 2*^0»5 cojttains some nultiterminal components then the matrix 
s 

may not he diagonal. If, hovevor, contains only two-terminal 

components then the ith component of E^ depends only on the ith 

s 

component of Y and therefore E^ is diagonal. 

2.11 Additional Remarks 

In sections 2.8 taid2.9 we have indicated the conditions 
under which the solution of a network can be ensured. These conditions 
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are mortJ general than those appearing in earlier work? this 
improvement is a conseq,uence of the discussion on the inclusion of 
dependent drivers given in those sections. Based on this discussion 
the development of a computer programme is currently under progress, 
which on successful completion is expected to have broader applica- 
bility than existing programme-s for similar analysis . Specifically, 
negative resistance characteristics are incorporated in the 
computer analysis xmder conditions vhere solutions are assured. 
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5. Sensitivity Studies, 

5.1 Introduction. 

In this chapter ve consider procedures for determining the 
additional information about the network regarding the changes in the 
response variables with respect to 

(a) component parameters, and 

(b) initial conditions. 

The models for sensitivity studies that are developed here 
admit very useful network interpretations. Consequently, the analysis 
programme that has been referred to in the previous chapter can also 
be used for sensitivity studies. 

Although time-domain notations are used throughout in the 
following study, it may be noted that this fact has nowhere been used 
for the development of the models. Consequently, the models developed 
are equally applicable to time- and frequency-domain analyses. The 
nature of the models is convenient for use in both analytical as v;ell 
as computer analyses. 

lirst-order sensitivity models were indicated by Leeds(9) • 
Similar, but more general models were derived for state equations (20) 

The formulation procedure given here is convenient for 
computer analysis of first-order sensitivity. The procedure has 
been extended to higher-order sensitivities . A simple recursive 
algorithm is possible for the linear case. 
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I'e shall consider time -invariant networks. Linear and 
nonlinear networks are discussed separately. 

5.2 Graph of the Sensitivity Model , 

In the previous chapter, we have discussed a procedure 
for determining the response of a network, which consists of the 
voltages Y(t) and currents l(t). Now we turn our attention to the 
variables in the- vectors 


(t) and (t), and also 

dP. Sp 


3 Y(t) 

o* P'1 ^ P 2 • • * S Pjc 


and 


i(t) 

3 Pi ^P2 ^Pk: 


In the above expressions, the symbols p, p^ , p 2 ,...,Pj^ denote 
the parameters of the network, dhey may also represent initial 
conditions of the nety/ork, In the linear case, they could be values 
of the components, i.e. E, L,. C, etc. The component parameters for 
nonlinear components will be defined later. 

For both linear and nonlinear networks, we have the f- 
cutset and f-circuit equations, 


and 


A I(t) =0 , 

B Y(t) = & 


3.2.1 


Differentiating Iq. 3 -2.1 parti ally with respect to a 
parameter p of the network, we can OTite 5 

.= 0 


3 p 


and 


B 


8 Y(t) 

3P 


= c 


5.2.2 
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The similarity of form between Eqs. 5*2.1 and 3*2.2 indicates 
that if components Virith appropriate characteristics are used, then 
it is possible to obtain the sensitivity coefficients 

as the solution to a network having the same graph as the given 
netx'-ork, 

5.5 Sensitivity Model for Linear Networks : 

We shall consider the sensitivity model for each type of 
component separately. In the end of -the section, we shall sum up 
the results in the form of a table . 

Independait drivers. 

The terminal equations of independent drivers are ■ 


= v(t) 


5.5.1 


for a voltage driver, and 
i = i(t) 

for a current driver. 

Prom these we get 

^(t) - 0 

a 


5.5.2 


and 


^ 1 

dP 


(t) = 0 . 


5.5.5 

5.5.4 


The counterparts of independent drivers in the saasitivity 
model have termiml equations of the form of Iqs. 5.5.5 an-<a 5.5*4. 
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In other words, the counterpart of a voltage driver is a voltage 
driver of zero value (a short circuit), and the counterpart of a 
current driver is a current driver of zero value (an open circuit) . 
Resistors. 

The typical resistor characteristic is modelled hy the 

equation 


v(t) = R i(t) 


3 - 3.5 


Differentiating Eq. 5* 3 .5 partially with respect to a 
parameter p ■?diich is differaat from R and G = 1/R, we get 


lZ(t) = R-^(t) 

dp ap 


3 .3.6 


Thus, if the parameter p does not refer to the resistor of 
-Q.. 5.5.5f then the ccunte apart of a resistor in the sensitivity model 
is a similar resistor. 

Next, let us find the counterpart of the resistor if the 
sensitivity terms 


dj(t) 

8r 


aid. 


SR 


are to he found with respect to the parameter R, Prom Eq. 3.3.5 we 
obtain 


§f(t) =R||(t) 5.5.T 

A nonlinear resistor can be described by a terminal equation 
such as Eq* 5 *5^7 • Thus, usir^ nonlinear network analysis methods, 
the resulting network can be solTed. 
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Alternatively, we can. revn:ite IDq, 5.3*7 a-s 



5.3.3 


I'iquatifm 3 * 5.8 can. be modelled by a parallel combination of a 
resistor of value H ard a current driver of value i(t)/R, T-ie 
direction of the driver is reverse of the direction in which the 
current i(t) is measured in the original network. The combination 
modelling liq. 3.3.8 is shown in Figure 5.5.'l(a')" 

For finding the counterpart of the resistor in the case 
where the sensitivity co efficients 


3G 


and 


3G 


are to be determined, we first re-write the terminal equation of 
the resistor as : 


i(t) = G ■v'(t) 


5.3.9 


Differentiating Eq. 3.3.9 ftlth respect to the parameter G 


v.'e get 

■, ■ ^ 3 ^ - ° 

As in the case of the nKodel discussed before, Eq. 3*5.8 
can be treated as the terminal equation of a nonlinear resistor. 

Alternatively, it can be modelled by the series combination 
of a voltage driver of value v(t)/G and ^ resistor of conductance G. 
The direction of the voltage driver in this model is the reverse 
of the direction in which the voltage v(t) is measured. The model 
is shown in Figure 3 . 3 . 1 (b). 



dR 



(b) 


Fig. 3.5.1 


Sensitivity models 


of a resi! 



Fig. 3.3.2 Sensitivity models of an inductor 


c+ QM 





Inductor 
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For an inductor we have 
di 


v(t) = L- 


dt 


3.3.11 


When the sensitivity pertains to a. parameter ■which refers to a 
component other ■than the inductor of Eq. 3 . 3*11 » we have 


_ L ^ ,Y 

3 P ■ dt ^ 3P ^ 


3.5.12 


In this case the counterpart of the inductor in the 
sensitivity model is a similar inductor. 

In addition to the terminal equation of the inductor 
in the sensitivity model, v^e must also determine its initial 
condition for a complete specification of the sensitivity model. 

Since p does not refer to the industor of Eq. 5* 3.11 > 

we have 


9i(o) 

3 P 


= 0 


3.3.15 


Hence, the initial condition on the counterpart of the 
inductor in the sensitivity model is zero. 

If the coefficiaits 

||(t) and |I(*) 

are to he determined, we have from Eq. 5.3.11 > 


5 v d / d i \ ^ 

■p “ ^ dt ^ a i ^ ^ 


dt 


5.5 ,1 4 


Equation 3 . 5.14 oan^ rewritten in the form 
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Equation 5. 3 *16 can "be modelled by a parallel combination of sn 
inductor and a cuourent driver, shown in Figure 3 . 3 . 2 (a). 

For the coefficients 

■^Ct) atd ||(t) 

(■where T = i/l), a similar analysis can be carried out to arrive at 
the model ^own in Figure 3 *3 •2(b) . 

Initial Conditions ; 

We must determine the initial -condition on the cohinterpcrt 
of an inductor in a sensitivity model in order that the latter model 
is completely deterministic . '.'Then the value of the inductor is 
varied infinitesimally, one may decide either to hold the initial 
current invariant, or to hold the initial flux invariant . The initial 
condition on the induG-tor in the sensitivity model depends on this 
decision. 

If the current is invariant, then clearly we have 

= 0 



B L 


3.3.17 




5«5*5 Sensitivity models of a capacitor. 



(c) 



Fig. 5 . 3.4 Sensitivity models of dependent drivers. 




for the model in Figure 5*5»2(a), or 
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_ 0 

B T . 

for the model in figure 3.5.2(h). 

If on the other hand the flux is invariant, then we 


5-5.18 


have 


Hence, 


^(0) = L i(o) = constant. 


3 i(o') _ Uo) 

31 - - L 

for the model of Figure 5-5-2(a). 

For the model of Figure 5.5.2(h), v/e have 

i(0) = T jZ5(0). 


5.5.19 


Hence 


^(0) = iM 

d T ' T 


3.5.20 


Capacitors. 


The typical terminal equation of a capacitor is 


i(t) = C 


dv 

dt 


5. 5.21- 


It can he shorm as in the case of components discussed 
earlier, that the counterpart of the capacitor of Eq. 3.5.21 in 
the sensitivity model, ■when the parameter does not pertain to the 
capacitor, is a similar capaci'fcor. 

When the coefficients 



are to "be determined, ive have 




dt 


f S-Z > 

^ ac 


dt 


= C 


L. 

dt 


B V 

a c 


+ V 


5.5.22 


Equation 5 • 3*22 can he modelled by a iqonlinjar capaoitor 
when rewritten suitably. Alternatively, it can be modelled by the 
series combination of a capacitor and a voltage driver, shov/n in 
Figure 5 ‘ 3 « 5 (a) . 

’//hen the coefficients 


(where S = l/c) are to be determined, the model shown in Figure 
5.5*3(’h) can be derived in a similar manner. 

As in the case of an inductor, the initial condition is 
required for a complete description of the sensitivity model. 

fihen the parameter p does not refer to the capacitor 
of Fq. 5.5.21, we have 


= 0. 3.3.23 


YJhen the parameter refers to the capacitor , two 
poasibilities aciso. Y/e may desire to keep either the initial 
voltage or the initial charge constant. If the initial voltage 


is invariant, we have 



-for the model of ^'‘igure 3.3.3(a), and 

^-YXQ). _ 0 

for the model of figure 3.3.3('b). 

’/’'ben the initial charge is invari;ait, ve have 

v(0) = = S (i(0) I 

aind hence 


v(0) v(0) 

9 'c " ” c 


for t’ne model of figure 3*3. 3(a), and 


3.3.25 


5.5.26 


d v(0) v(0) 

SS " S 


for the model of figure 3.3.5('b).. 

Consider the case vj-here the initial current through an 
inductor or the initial voltfge across a capjscitor is the parameter. 
In f/.io case, -j^e can derive the terminal equations of the counterpart 
in'.luctors and capacitors from iq's. 3.5.11 and 5.3.21 respectively. 
iL'hese counterparts are linear inductors and capacitors identical 
to the ones in Eqs. 3.3.11 aid 5.3.21. The initial conditions on 
these counterparts are a unit current and a unit voltage, respective! 
since 



3.3.27 


for an inductor, ".nd 
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S y(0) 
a v(6) 


3.3.28 


for ot capacitor. 

Dependent drivers. 

We restate the equaticms of the four types of dependent 
drivers here for convenience. They are ; 


Type 1 

VgCt) = v^(t) 

T’rpe 2 

V 2 (t) = ^ i^(t) 

•Type 5 

i2('t) = "y 

Type 4 


3.3.29 


5.5.30 


5.5.32 


■'.'-'hen sensitivity coefficients o.ith respect to a parameter 
not referring to a dependent driver are to 'be determined, it can "be 
sh't.n that the counterpicrt of a dependent driver is a similer 
d.epend.ent driven ■'■'e shall discuss the models of the counterparts 
of dependent drivers in jensitivity models for their om parameters. 
Type 1 


From Fq. 3.3.29 iwe have 




= c< 


Hi 


(t) 




3.5.53 


■quati n 5.3.33 can 'be represented as the series combinaiicn of a 


dependent driver and an independent driver of value v^ . Such a model 
is shewn in Figure 3.5.4(a). 
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Sensitivity models for dependent driver types 2,3 and 4 can 
be derived in a similar manner. These models are shovn in Pigires 5.5 .4(b 
c,d) respectively. 

ihe results of this section can be summarized by the following 
algorithm ; 

Algorithm 5* 5.1 

lor obtaining a netvork "vshose rest'onse consists of the 
sensitivity coefficients 


B 7 
B p 


(t) and 


dP 


(t) 


yhsre p could assume any one of the -arameters R,L.C,T,S, G, o4 , , If , 

or S , or one of the initial conditions i(0) or v(0) , go through 
the folloTA’ing steps s 

(i) 'Short-circuit all independent voltage drivers and open- 
circuit all independent current drivers of the original network 

(ii) Connect the appropriate driver in parallel or series 
with the component tdiose parameter is p, and / or set the appropriate 
initial condition on that component as indicated by Table 5.5.1. 

Set 0.11 other initial conditicns to zero. 

5.4 Higher-order Sensitivities. 

General higher-order saasitivity coefficients of the type 


s Pi a P2— 3 It a Pi Pg - - a Pfc 

v/ill be discussed in this section. On the basis of the results 
obtained in the last section, we can straightaway conclude that 
independent voltage and current drivers in the original network axe 
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First-oxder SensitiTity models for linear components, 


Component 

Driver 

Initial Condition 

Type 

parameter 

Type 

Placi ng 

Value 

i or V q or 

invariant invari ait 




Capacitor 



Dependent 
Driver 
Type 1 


T;;^pe 2 


Type 5 


Type 4 


Current Parallel i/E 


Voltage Series v/g 


Current Parallel i/L 


T 

Voltage j Series v/T 




Voltage Series v/c 


Current Parallel i/s 


Voltage Series v 


Voltage Series i 


Current Parallel v. 


Current 






-io/L 


io/i 




-vo/C 


vo/S 











































replaced "by short :ind open circuits, respectively, in the sensitivity 
net..jrk. Similarly it is quite ohviais that if the parameters 

- - _ - do not refer to a certain compjonent of the R, C, or 
dependent driver types, then the counterpart of such a component 
is a. similar component 5 and in the case of a capacitor or an 
inductor the initial condition is zero. 

Te shall now derive models for the cases when one of the 
p.'.ramete :s refers to the component. For example, in the case of a 
ro sis tor whose terminal equation is 


v(t) = R i(t) 


3.4.1 


Sup"ose it is desired to find the sensitivities 




^ a P. 9 p, 


and 




aR i)P2 - — 3 p^ 


where none ■'f V2>~ " ~ to the resistor of hq, 5.4 .I. '^e 

have from Fq. 5 » 4.1 " 




R 


c* 




BR^Pg 0K 8^2 


eJPo 


•C* Pi, 


5.4.2 

As a more general, case, consider the following sensitivity coefficients 
^ Y B I 

— — h.; : — — — ^ — — — — ^ — — p 

^ ^ ' ^ ^ c;- ^ ^ Pkl+I ^ Pk 

in which the parameter R appears k- times . 
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Again, from Fq, 3.4.1 "we hs;ve 



•g) i . 
*dP2 


+k. 


k-i 

B i 

S R’‘aP2-BPk 


5.4.5 

Equations 5.4.2 are a special case of 4q. 5.4.5 where = 1 . 

Equations 5.4.2 and 5*4.5 are similar to Eq. 5.5«7 af 
the previcnis section. The model shown in Figure 3.5.1(a) can be 
\j.sed to model £qs. 5*4.2 and 3.4. 5» when appropriate values of the 
driver are used. Similar equations and models can easily be 
dsrived for Inductors, capacitors, and depaident drivers. Models 
slvown in Figures 3*5.1 to 5*3.4 can be used with the values of 
•drivers given at the end of this section in Table 5*4.1. 

Initial conditions ; 

Here once again we have the choice of holding either 
the flux or the current invariant (in the case of an inductor). 

If the current is held invariant then the initial 
cond.iticn on that inductor is zero. Similarly, if the voltage of a 
capacitor is held invariant, the initial condition on fiat capacitor 


is zero. 

If the flux of an inductor is held invariant, and 

if sensitivity coefficients of the kind 


T(t) 


cond 


a^i(t ) 


are sought, then we have 
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^ \(oy )^o (-0^ k . i(o) 


If the sensitivity coefficients are sought vAth respect to any 
parameter other than L, the initial condition on the inductor is 
zero. Similarly, if sensitivity coefficients of the type 



and 


ac’" 



are sought, and the charge of the capacitor is held invariant, '.ve 
Iv’-ve 


B \(o) 



( f (0) ) 


(-1)^ k . v(0) 


5.4.5 


If sensitivity coefficients with respect to T in the case 
of an inductor or S in the case of a capacitor are sought, we have 
a simpler situation, v/e have in particular, 


^ ^ i(0) 


0 f or k > 1 ; 


and similarly 


.0 f or k > 1 . 


In connection mth sensitivity coefficients of the t-^rpe 




^ P l - - - 


and 




B Pi ■ " “ 


3 . 4 .6 


5.4 .7 
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k if one of th3 parameters is an initial 

condition, the corresponding inductor or capacitor has no driver 
associated with it and a zero initial condition. In other fiords, 

Jhe sensitivity model of that componait is just like tie case where 
tie pp-rameters do not refer to that component. If two of the 
paraneters are initial conditions, then the coefficients sought 
are all zero. 

In the above analysis, we have found a method of finding 
tie kth order sensitivity coefficiaits based on the knowledge of some 
coefficients rf order (k-l). However, this does not pose a serious 
problem since the first-order coefficients can be found with the 
knov, '■ledge of only the network response, 'and then the second -order 
coeffici ent s can be found using the first-order ones, and so on, 
upto any desired order. 

Let us summarize this section again in the form of a table 
and on algorithm for finding 


a Pi &P2 SPi 


^ T 

and — where k 1 . 

^ P2 - 


Al.gorit’''m 5^4.1 (i) Short-circuit all independent voltage drivers, 

and. open-circuit all independent current drivers. 

(ii) Set the initial conditions on all dynamic components 
except the ones -whose parameters are in p^ , Pg, ...jP^ ,to zero. 

(iii) Connect the appropriate drivers in series or parallel 

with the compon-nts whose parameters are in p.^ , . . . .,Pj^, and/or set 
the appropriate initial condition on that component, as indicated 
by Table 5.4*4* . 
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iabli- 5*4»'I Higher-order Sensitivity models for linear components. 


Compcnent 


Associated Driver 




Initial Condition 



T j Voltage | Series 


k ^k-1 
1 d V 

T V3P2**-aPk 


0 (-1 ) (k-l) 


I 

i 0 if'k>1 



Capacitor 1 C | Voltage | Series 



_ I (-l)^(k-l),o 


:D apend ant 
Driver 
Tvtdo 1 


jType 2 


a P2 • • • aPk I 


S Current j Parallel — r -rr — i 0 if ,k=1 

S SPv I S 

, ^ 0 if k > 1 


Voltage j Series 


Voltage j Series 


k. .k-1 

_L a ^1 
3P2'*-3Pk 

, ^ i 

> *" , 

1 1 

! 1 

i ■ ' 

- 

k ^k-1 . 

! . i 

i ^ 




jType 4 


Current Pars 




h a’"-' i, 


3P2*‘^* 

* 

1 




















5-1 


example 3 *4 -I 


For the triode amplifier shown in Figure 3*4 3- find 


First, 


A - ^ 4 d A 5^A 

■'i ’ -au ’ 3\ > 3 b 2 • 

L 

to simplify, tfe connect a unit voltage driver at the 


in- ut . 
i’hen, 


V 


2 


A. 


The equivalent network is shown in Figure 3.4.11. The sensitivity 

yi jA 

models for = , and — -r are shown in Figure 3 4*1 c,d, and e 

BR dw ^ 

L 

respectively . Solving these we obtain 

A = - U 5 




3 A 

3^ 




R + R_ 
P L 


d A 


E 


- 






and 


2 


V * V' 


3.5 Ronlinear Network. 

In the case of line.ar networks, the component values were 
used as parameters of the network. No such unique and natural 
association is possible in the nonlineox case. For nonlinear networks, 
a parameter (always denoted by p) is introduced in the terminal 




Fig. 5*4 •'1 • Sensitivity of a triode amplifier 
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equation of each type of component other than independent drivers. 
Pile relation of this parameter to the function in the terminal 
equatirn 5. s not prescribed. The sensitivity coefficients are then 
doscribed in terms of the parameter p. Consequently, the models 
involve the partial derivatives of the functions in the terminal 
aquations of nonlinear components ivith respect to p . The 
component characteristics ajce modelled by the following equations 
when tliD parameter p is introduced ; 


liesistor ; 

i = G(v,p) 

3.5.1 

Inductor ; 

d^ 

^ " dt 

3.5.2 


i = T(0,p) 


Capacitor 

i S3 - 

dt 

3.5.3 


V = S(q, p) 


Dependent 

Dr iy er s 


Type 1 

V2 = V (v^,p) 

5.5.4 

Type 2 s 

Vg = V(i^, p) 

5.5.5 

Type 3 . 

i2 = I(v^, p) 

5.5.6 

Type 4 r 

ig = l(i|» P) 

3.5J 


Prom these \V8 obtain the sensitivity models as follows : 



Resistor s 
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5 i 3 Q a V 
cvp' ~ 5^^* c5 p' 


3.5.8 


vrhere p' does not refer to th't resistor. 

For its own parameter, we have 

^ i ^ G '5 T 0 Q 

B P ” * 9P B-P 

In the atove models we have seen that the knowledge of the 
derivrttives 


5G . ^ G 
'-r — and — 
dr dp 


pertaining to the characteristic of the resistor are required for 
its sensitivity models* If second -^order sensitivity coefficients 
■were sought, we have the following equation s 

0 ^ g ^ ^ Ll Ijl di. _iLL- ^ 

a P^P' “ 0^ p, ^ * ^P ■ Bp' "^dPa. ^ BP ’ 

5.5.10 

Tims, for each increase in the order of the sensitivity coefficients, 
additional information about the component terminal equations is 
required for the sensitivity model. 

T'le shall abandon the study of higher-order sensitivities 
of nonline<tr networks, since it is impractical to require so much 
information about any network vdiich is big enough to justify computer 
solution. However, there is no theoretical difficulty in constructing 
models for higher-order sensitivity coefficients. 
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Let us continue the discussion on first*-order models, 

Por an inductor, ?jhen p* is a parameter not refer ing to it, 

■'■"e have 



f ^ ) 

dt ^ 3p' ^ 


^ i(t) 

p' 



Mil 

d p‘ 


And for its own parameter p, 


and 


'Ll. . ( 11 

d P dt ^ 9 p 


) 



3.5.11 


5.5.12 


For a capacitor, for p' which is not its own parameter, 


' 5 > i d / o> q N 
/9P' == dt A 9 p' ^ 

and ' 3.5.15 

5 V 0 S 9 q 

9p' O'! * 3p' 

And for its own parameter ; 

■9 p ' dt > 9P ^ 

and 3.5 .14 

9 V O S 3 S o* q : 

Op “ ap ^ 9<1 * 9 p 

For the fcur t^rpes 0 f dependent driver s we obtain the follovsdng 
equations, "Aiere p is the parameter of the driver cone emsd and p' 
a different parameter. We have t 
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Typs 1 ; 

■^p' ” ap' 



^ Y 

9^1 

or ^ p ‘ 

“Sp 

* 9P 

Type 2 ; 



9 ^2 

a V ”3^1 


3p' ' 

a p' 



By Sy 

ill 

or ^ = 

o)P 

Op Oi, ' 

ap 

Type 3 ‘ 



3^2 

3 1 a''! 


'ap ' ” 

* 9p' 


3 "2 

3i 3i 


or = 

^p 

9 P ‘9'^i 


Typ e 4 i 



3 ^2 

^ "2 


II 

gi^ * ^p' 


*5^2 

9 1 Si i 

9^2 

or rr — - 

ap 

c9p ^ aia* ' 

SP 


Initial Conditions : 


Indnctor ■ If the flax' is held invariant , then ve have 


.3-5.15 


3.5.16 


5.5.17 


5.5.I8 


5.5.19 


5.5.20 


3.5.21 


3.5.22 


aP 


5.5.25 


If we try to hold the current invariant vre get 
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or 


^ i 

') T "3 T 


■^p“ 


a p 

p 

^ 0 



^ p 




3. .5. .24 


3 . 5 .25 


. T / 

provided —.^ ^ Oj.ishich nay not be true in general,. 

If is treated as a parameter and the sensitivity coefficients 
with respect to it. sought, then we have 


'2 P 


3.5.26 


ind the terminal equation is valid for that inductor. 


Capacitor ? 

As in case of the inductor, it is not possible in general 
to hold the voltage invariant. If the charge is held invariant, 
"■re have 


2% 

9P 


5.5.27 


and provided dS/dq^/O to can hold the voltage invariant by making 


d V 


9 S 


'bS 2% 


■ 3 P D<io* P 


= 0 


3 .5.28 


■9 '1. 


or , 


K) ■ 


-9 S/'^ p 

■'Wit 


3.5.29 


If q is treated as the parameter with respect to which the 
^o 

sensitivity is sought, then 


9 q 


( 0 ) = 1 


5.5.50 
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and the terminal equation 3»5»12 holds for that capacitor. 

In conolusion, first-order sensitiTity coefficients for fixed, 
nonlinef.r networks can he ohtaiix^d as the solution to a time-varying 
netvsrk, all hut one of whoso components are linear, and the 
roiaaininr one has a straight-line terminal characteristics not passing 
thivjiigh the origin. 

The following algorithm is a collection of the reailts of this 

section for purposes of apT.Ii cation. It results in a network 
, c)iV ^I 

v.'aoBu solution is and ^ , vihero' p is the parameter of 5qs. 3*5*’' *• 
3.5 «7 or rxL initial flux or an initial chnige. 

Algoritlim 3* 5*1 (i) Short-circuit all independent voltage drivers 
and open-circuit all indepen^.ent current-drivers.- 

(ii) Replace all components except the one whose parameter 
is p hy the one indicat sd in Table 3«5»1 ond set initial conditions 
to zero. 

(iii) Replace the component whose parameter is p and adjust 
its initial condition by the one indicated in Table 3* 5-2. 

.'ilxcmple 3 . 5 •'I 

It is desired to study the rate of change of the amplitude 
of oscillation in the net ,..ork of Figure 3.5-1 a with respect to 
a change- in R. The choract eristics of the tunnel diode G is given 

as, .'1: :■ ' ' ■ ■ 

i = -.02v+ 0.00002 v^, ;3*5*31 

where i is in milliamperes and v in millivolts. 






Table 3.5.1 
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1 Component 

L- • _ 

Linear 

Replarenient 

! 

i Type 

} 

Terminal 

Equ ati on 

No. 

Type 

Terminal 

Equation 

No. 

j 

1 'Resistor 

5.5.1 

, 

Resistor 

■ 

5.5.8 

jlninctor 

1 

3.5.2 

Inductor 

5.5.10 

1 

'Ca 'acitor 

! 

5.5.5 

^ ^ r 

! ^ ■ j 

Capacitor | 3*5.12 

I i 

jDepenicint 
'Driver 
jT.ypie 1 

5.5.4 

i 

Dependent 

Driver 

Type 1 

^ ^ 

5.5.14 ! 

I) opend 3nt 
Drivar 

Type 2 ■ 

1 5.5.5 

Dependent 

Driver 

Type 2 

5.5.16 i 

Dependent 

Driver 

T'^8 3 

3.5.6 

Dependent 

Driver 

T'.'pe 3 

5. 5.18 

Dependent 

Driver 

Type 4 

5.5.7 

■! 

Dependent 

Driver 

Type 4 

■ 

3.5.20 

1 









Table 3.5.2 
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Compcn^-nt 

Replacement 

Initial Condition 

i 

jTvna 

I 

I 

T.D. 

Do. 

F. am- 
meter 

i 

Type 

1 

T.E. 

Do '. 

Q or 0 
^0 “0 

const^mt 

i or 

0 0 

consta-- 

nt. 

} q or 0 
! ^0 ^0 

1 para- 
j meter 

i 

Resistor 

3.5.1 

i j 

p 1 Resistor ) 5 . 5.9 

i j 

- 


- 

Inductor 

} 

3.5.2 

— 

P 

: 

Inductor 1 5 - 5*11 

i 

I 

0 

' 1 

^ o7 

0 


5! 

n 

^0 

n 

3 . 5 . 1 * 

1 1 , 

'! ' I 1 

Capacitor 

j 

5.5.3 j p 

1 

1 

C apac itor 

5.5.15 . 

. 

0 

'bs/y 1 . ! 

TsTTq” I i 

^ ^ 1 ^ 1 

tt 

ti 

a 

ft 

3.5.12 

{ " ' i 

-1 - 

Dependent 

Driver 

Type 

! 

3.5.4 j P 

1 

Dependent 

Driver 

Type 1 

3.5.15 

- 

- 

- 

Dependent 
Driver 
Type 2 

5.5.5 

p 

Dependent 

Driver 

Type 2 

3.5.17 

■ 

1 

- 

- 

- 

Depaidenti 
Driver 
Type 5 

* 

i 

3.5.6 

p 

Dependent 

Driver 

Type 5 

3.5.19 



■- - ' , 

Dependent 
Driver 
Tyne 4 

5 .5.7 

■ 

p' 

Dependent 

Driver 

Type 4 

5.5.21 
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I'igiire 5 *5 •’1 ("b) shov?s the sensitivity net 7 ;ork ivhose 

solution is , 

(iR 

In Figure 3.5.l(-b) , 

G-^ = - .02 + .00006v^ 


5.5.32 


■e ohsarvo that at the instant 


t when 

Jr 


"'peak. 


rrj~ = 0 and therefore -rrr 

c) ^ ^ it 


~S\eak 

SR 


Thus, siffluntaneous solution of the not’®rks of Figures 3.-5.ta and h 
will yield the desired result, 

3.6 i'loasures of Sensitivity, 


■'fe have so far considered models for sensitivity analysis. 
Sensitivity analysis constitutes an important link hetvireen the analysis 
■and design of networks. In particul.3r, with the advent of integrated 
circuit techniques, design ohjectives "based on sensitivity have hec»me 
valuable ( 14 ) • 

For basing design objectives on sensitivity one needs non- 
negative real numbers based on sensitivities, called measures of 
sensitivity. In this section, measures of sensitivity based on both 


time-as woll as frequency-domain sensitivity coefficients are discussed 
A measure of a single coefficient; 

Consider the real number s defined by : 
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where u is a variable (voltage or current) of the network and p is 
a pirametor . This number is an upper bound on the relative change 
in u(t), for any time t, for a small change in p. In other words, 
far small charges in p, sSp is a close approximation to the 
change S u in u, whore % is the chEinge in p. 

\fcje can easily define a similar quantity, a vector S, for a 
vector U of network variable*, and a single parameter p as s 


S 


max I ^^1 ! 
I ci P i 


i 

max c* n j 

t : 


3.6 .2 


1 voightod scalar moasare carrying the information about m parameters 

,is 


m 

j = w. s. 5 *6 .3 

7^ 11 
1=1 

where are the weights associated viith the parameters p^ and 


s 


for a scalar variable u. 

If the parameters p^ are knoxm to vary only within p^ + ki > where 


= max 
t 


j 'B u 

I dPi 


5.6.4 


then 


A i 


<< 


Pi, 


■ ■ xmT 

X^ermissihlf; 

chan^es 


C u 





wherd • Thus, the measure J-in E< 1 . 3 *6 .5 is a measure of 

the worst-caee change in u for restricted changes in 
The measure 


K - w. s/ 

1=1 1 i 


3 . 6.5 


can similarly he given an interpretation. 

If p^ have a Gaussian distribution with P- > then we 


<- d / 0\x ^ 


2 / "b u \2 
^ BPi ^ 


3 . 6 .6 


Hence, if w. = a t then we have 

X X 


6 K. 

^u '' 


3 . 6.7 


In In words, K is an approximate upper bound on the variance of the 

2 

response variable u, vhen the variances of the parameters p^, 

are smell and knoon. 

t’requcncy-Pomain Measures for Linear Networks : 


Consider a variable of tte network in the frequency 


domain 


U(s) = 


gi - h) 

,4 ■ b) 


3.6.B 


The locations of polos and zeros, that is r^, . . . . ,r^ and z^ , 
are the important factors govering Il(s) . The folio -wing real 
quantities one meaningful measures of frequency-domain sensitivity 


1 


5 . 6 .9 
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a 


3 p 


3.6 .10 


m 


j=1 


ill 


hk 

3 P 


n 

j=1 






9 p 


5 . 6.11 


C’lj third measure is actually a weighted sum of the first two, with 
the (different components getting equa,l vK.ightage on the loga,rithmic 
scale of frequency. 

'fho above measures are for a single parameter p. They can 


he extended to several parameters either hy the worst case method 
or by "tte statistical method (See Sqs. 3»6.3, and 3*6. 5)* 

3.7 Stability analysis by sensitivity methods. 

In the Liapunov sense, the motion of e. system described by 

I 

o 

X = p(x, u) 


x ( o ) = x ^ 


3-7.1 


is called locally stable if for any > 0 there exists a. 


0 


such that 


X - X’ | < -■=^ I X(t) - X'(t) < f V t 

o 0 ) -^ t I 

h^ere X' (t) is a motion of the same system, governed by 


X = P(X, U) 

starting from the initial state 
X(0) = X»^ 

(see (17))- 


5 . 7.2 


3.7.3 


5.7.4 
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If tixfe coofficicuti 


s . . = mx 
10 ^ 




d* 


3.7.5 


ara all finite, then for 


S. - 


€ 


^ ^ s. . 

1 0 10 


3.7.6 


the con'iition 5*7 *2 is satisfied . If iiny s. . is infinite, then the 

1 J 

00 ni it ion fails for the initial state 


X' = X^ + 
o 0 


0 

0 

€! 

0 


5.7.7 


where ‘6 is in the jth place and is arbitrarily small hut non-zero. 

Thus wa have proved that the cond-ition that s. . are all finite is 

13 

equivnlont to the condition For finite inteirvals of time, s^^ 

and their combinations indicated in the previous section sem/s as 
moasur s of instability of the system. 

5.8 Further 'Rem.erks. 

Sensitivity coeffiaients of an arbitr;:ary order can be obtained 
for a linear nettaork, given only that amount of information which is 
nerded for its analysis. For a nonlineto network, however, additional 
informs-tj-on is required for. each order of differentiation. A 
•DroceAu;ce is established only for first order sensitivities in the 
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case of nonlineiX net'.;orks, for it seems neither practical nor worth- 
whiloi to ask for the laige amount of infomation ’Aiich. is needed 
for higher order coefficients. Thera is, ho\vevor, no theoretical 
difficulty in arriving at these. 

It seems appropriate at this place to note that Leeds 
and Ugron (10) have indicated a method of finding the sensitivity of 
a single variable of the network to all its parant^ters by solving 

only ono sensitivity network. The method is useful for only linear 
net'TOrks. 

The question of measures of sensitivity \d.ll be further 
discussed in the next chapter in the light of a design problem. 
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4. Opt isdz at ion Studies, 

4.1 Introduction. ‘ 

\iJe have so far developed models for sensitivity studies, 

essentially another facet of analysis isiiich extends the range of 

possihlfi use of the existing methods of analysis. Our concern nov» 

is to direct these analytical studies for design ^purposes. In this 

context, optimization studies constitute an essential link hetween 

analysis and design. For example, one might he expected to select 

a host network from a continuous class of similar networks based on 

an ap.'Topriate design criterion. That is, a propsr choice of 

■naramjtors of the network within the specified hounds satisfying an 

ohj octive critericn is an important facet of design studies. 

Optimization studies, such as the prohlan of determining a 

network satisfying the conditictn of sensitivity minimization have 

received attention in the literature (1,14). The very nature of 

these design studies is such that they do not suggest unique methods 

t 

of solution. 

In the following section we shall formulate the parameter 
optimization prohlan for a linear netv;ork and a general design criterion 
The formulation is in terms of the values of the components ; it leads 
to the soluticn of the prohlem as a succession of linear programming, 
prohl ons . Some adviantages of the use of component parameters as the 

basis for optimization studies are pointed out. , 

Another prohlem that is discussed is that of optimal control 
of a network. In other words, optimal values of the indepei^ent 
drivers in a given net?/ork are determined as functions of time on the 
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basis of .an appropriate design criterion. This problem has been 
thoroughly discussed in control theory, but without reference to the 
int'-rnal structure of the system (i1,18). 

'J^'he solution to such a problem can be obtained, in terms 
'.f .a t’WD-point boundary value problem for .the state equations of the 
network, .-ind another set of first-order differential equations of 
similar dimension derived from the state equations, ’which .are called 
co-state equations of the network. 

W'e show that the co-state equations admit a simple network 
interpretation, which is useful in systematizing the formulation of 
these equations. Besides, it may also be possible to use network 
analysis programmes such as the one referred to in chapter 2 for the 
solution of the resulting two-point hound .ary value problem in some 
cases. These, possibilities are explicitly pointed out. The theoretical 
bac’^ground required for the development of the model is discussed in 
a separate section. 

In some cases, the solution to the problem h.as to be further 
testadl for optimality. Sufficiency conditions for optiimlity have been 
discussed by Robbins (21 ), It is demonstrated that these conditions 
do not depend on the gro-rh of the network and therefore do not admit 
a net T^ork interpret .,i.t ion. 

4.2 Parameter Optimization - A Component -level Formulation. 

of network synthesis do not t.3.kB into 
account the effect of non-ideality of components on network behaviour. 
The main emphasis in synthesis procedures has been to realize a network 

having specified nominal performance with as few components as possible. 
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In large networks, differences be tureen individual component 
values nay have ccnsidarable effect on the response which may become 
unacccptahle merely because of this reason , The problem is p.articularly 
serious in integrrted circuits, component replacement is not 

■oossible. On the other hand, the importance of the number of 
components is not great in integrated circuits. 

As a consequence of such developments in the technology 
tsf network construction, considerable interest has been shown in 
recent woxiq in associating sensitivity studies v/ith network synthesis 
(1,8,10,14) . A parallel dsveloj'msnt has also taken place in the 
literature on autom(xfcic control; .an extensive review of such development 
can be feund in ( 6 ). 

Schoeffler (I4) has emphasized 'the importance' of sensitivity 
■analytsis in the design of net'v/orics and has presented a method of 
sensitivity minimization. Bolstein (1) has formally analyzed 
ichoeffl'u: ' s approOLCh and indicated a more general class of problems 
of -which it is a special case. He has also suggested other criteria 
for optimiz.a.ti fn based on sensitivity. 

In this section, ve sh.all formulate the problem of the 
minimizatien of a general objective function. The formulation uses 
finit.3-difference equations describing the manifold in the parameter 
sp.ace 7.hich constitutes a set of continuously equival'ent networks . The 
■orocedure uses the values of the com-ponents as the variables for 
forTiulflirtion. 'rhis procedure is helpful in roducing the algebraic 
m.oniipul.ation to a minimum both in the formulation as well as solution 
of the problem. 
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This secb.on is concluded with an example to illustrate, the 
m-tiod. jJurthcrmore , the example slso indicates th.at it is spELetimes 
■oossihlo to improve the insensitivity of a netyork without increasing 
the nuniher of components; moreover it may he possible for a network 
containing p smaller number of components than the general member of 
a continuously equivrlcnt class to be the most insensitive one . Only 
line .-x net', orks are considered, 
formulation 

”v7e start from a network ''.“hich realizes the desired netwoik 
function md refer to it as the ■'sta.rting network". 

How consider all networks which have the same configuration 
as the storting network, o.nd \laich realize the s^ime net .-^ork function, 
but differ from the st-arting net'work, and from ejicb-other, in the values 
of their com’ onents, Ea.ch netnoik of this set can be represented by 
the vector of its component values ; 

P = col. (p^, ,p^) 4.2.1 

The collection T of all such n-vec tors, 

P = P I T(P) = ^ 4.2.2 

where T(P) is the specified network function, and its desired value, 
forms a mnifold in The above set-theoretic notation is used 
to confoma to the sto.ndard notation in tte literature. 

Por P to be non -trivial (i.e. for it to contain more than 
just one point ) the startir^ netwrk must have some built-in 
red-undancy. CFte shall return to this point when discussing the example.) 



Consider t-n neigMbouring points, P and P + SP, I- 
«tainir« only the first-order tern of the Taylor's expansion for 

tIp + 6p) 5 


T(p + S p) = N'H ’ ^ ^ ^ 


^.2.5 


where 


^ ^ ’ SPn' 


axid ^ / 


V denotes the scalar product, ttet is 


<x,y> 


T 


Since 


T(P + & P) - -*c 


4- • 2 •4' 


4 *2 -5 
4 *2 *6 


’e have, 


<:||.sp> -0 


4.2 .7 


ij-^g left-hand member 


of hq .4.2.7 is a 


function 0 f the Laplace 


this eonation vail give risi 


se to several 


vaxi able s, and hence 

constraints of the t,ype •• 

- 0 i = '* »2, •• • 

constraints -depends on 


</Pi(P), &P> 


4 . 2.8 

the complexity of 


Tbe number h of tbsse 

,(,). one «y of 
° X ' 


of constraints is to nrite 


4.2.9 



where iii^(s) and n^(s) are polynomials in s, having no common 
factors, and fbe coefficients of m^(s) are functions of p and 6 pj 
and then equate each power of s in m^(s) separately to zero. 

Another method of obtaining a set of constraints 


is to vrrite 

q 

a + a. s + + a s 

° ^ % 

T(s) = 4.2.10 

h^ + h^ s + .....+ h^ s'" 

where a = 1 

q 


and then write 


< 


3 \ 
”3 P 


, —0^1 — 0, 1, ...., q“’1, 


4.2.11 


and 


/■ 


\ 


\ 


3P 


, S P > = 0 *, i = 0, 1 


■» r 


4.2,12 


In arriving at Eq. 4.2.8 we have truncated the Taylor’s 
expansion for T(P + ^ P) after the first-order term. Hence, we 


must impose bounds on | 
reouireroents, typically s 


in order to meet the accuracy 


I ^ Pp I j^ ,i^1,... .^n 4 . 2 13 

The values of in Eq. 4.2.15 are based on a suitable 

criterion for the control of par-step error. 

For examv-le, we may require the design function T(s) to 
be held v/ithin + per unit in the useful interval of frequenoies 
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on the jw axis. In this case, the following mle can. he eaplojred ; 


A,. 


n X max 






bT( j») 




B P, 


i 


4.2.14 


for each i = 1,2, . . . ,n. 

Alternatively, we may decide to control the error in 
each term of the expansion of T(s) given hy Sq. 4-2.10 so that 
a certain accuracy is maintained at any frequency. In such a 
case, the following rule can he used : 


n X max 


1 


% * ^^i 






a Pi 


(-L .i^oj 

1 : 

|_L l!l 

ll 

h ’ t^p. 

0 1 ^ 

1 

ti ■ Sp't 1 

i 


In addition, we must enforce 

- Pi 


4 * 2*15 


4 .2 .16 


in order to maintain passivity. 

can comtine Eq^S'. 4.2.15 and 4-2.16 into one simple set of 
constraints *by means of the transformation i 


A . - min 


{pi> } 

rr . 

Sp* i= S Pi +4 i ’ 


4.2 .17 

4. 2. 18 


/ n 

Aj “ -Aj + A- 


and 


4.2.19 
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In terms 


Tis of tie new varieties S P' i 


4.2.20 


o4- ’ i= “I* • •'»”• 

In terms of the transformed variatles Bq.s. 4-2.8 are as folio -.s 

<F. (I). Sp'> = <ri(p). /t'' > 


4.2.21 


where 


gp» = Col. ( Sp’n ) ’ 


and 


. ( , -A 2> •• • •» An ' 

e to he minimized *. 

4.2.22 

minimize at each 


A = Col. V A-j » 2.\ 2 

I et he the measare of performanc 

0=0 (P) 

3teepest-4.c.nt t.chm<iue « can 


Then, using the 
step the change : 


or equivalehtlv 


S js = ^ ^ ^ ’ 


4.2.23 


U ' .■ 2JL , I: P' > 

At some, stage we shall observe 

min (S0) = 0 


4.2 .24 


or, 


min 


■v4iich ensures 


■ /M \ N 
(s 0‘) - \ ^ 2 

that a l>aan reached. 


4.2.25 
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At this point, a discussicn of the merits of the above 
formulation seems to be appropriate « sTe have throughout used the 
values of the components arranged in the vector P as the variables 
for the formulation. This approach eliminafes the additional 
algebraic manipilation- involved in the transformation of the 
results to the values of componaits, ishich will be necessary if a 
transformed set of variables is used for the formulation. 

An examination of ISqs/ 4,2.14, 4 » 2.15 will demonstrate 
tbot any error criterion can be directly translated into bounds 
on the step size of the variation of parameters for the purposes 
rf optimization. 

The forraulati® of the sensitivity minimization problem 
given in (I4) is based on the assumption that the net's'/ork contains 
a component of each type between every pair of nodes. Por 
example, in this treatment, in an LC network, the general member 
of the continuously equivalent class used for sensitivity 
minimization would contain ’an L as well as a C between every pair 
of nodes. The necessity for such elaborate configuration does not 
arise in the present formulation. 

Furthermore, the trajectory on P- chosen in (14) niay not 

contain the minimum number af constraints, and thus may lead to 
restricting the choice to a subset of P. In the present formulation, 
however, the constraints chosen are a necessary set of coiiditions 
for the invariance of T( 8 ) £ps.d heiffie are , nonrest riot ive u'ithin P. 

Newcomb ( 22 ) has raised the question as to whether 
one could obtain all possible realizati ons of the network function 
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’■•hich have the same graph as the starting network hy the process of 
choosing neighhoaring points on P. Since we have fornmlated finite- 
difference constraints which describe P in Eq. 4*2. 7» this question 
toils dovn to the question of connectedness of P. The manifold P 
nay he simply connected, like a circle, or may consist of more than 
one part, like a hyperbola. If P consists of more than one part, 
then 'q. 4«2,7 -111 clearly restrict the solution to. that simply 
connected part of P to which tte starting network belongs. 

Another p^int that must be raised in this respect is that 
of the exist aace of local minima of the objective function. If such 
minima exist, the desired absolute minimum may not be obtained 
directly by the procedure discussed above. The study of these 
properties can be made only for particular types of networks and 
objective functions. 

■'5'XBmple 4*2.1 

Consider the design of a realization of the driving- 
point impedance 



s + 1 s + 5 


4 . 2 .26 


A realization of Z(s) is given in fig. 4-2 .1 (a) . 

Thi s rcalizati'm does not eontain any redundancy, ’’e introduce a 
zero resistance R2 to obtain a non-trivial continuously equivalent 
class P, (See fig. 4 -2 .l(h) ) I'he network in fig. 4.2.1(b) is a 

special case of the non-series-parallel realization described by 
Lee (S)* In terms of the component values, or parameters, the 
driving— point impedance Z(s) is expressed as . 
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(a) 

A Canonical Realization of Z(s). 



The Transformed Reiv/ork. 

I’ig. 4.2.1 Sensitivity Optimization of an R-C network. 



s(p^+ P2) S2 + (Pj + P5) 


^ ^ p ^ + S P ^( p ^ + P ^) + P 2 ( P 4 + ^ 5 ^ ^ ^ 


4 . 2.27 


■ 7 , “He re 


3 ind 




F = Col. <C^, Cg, R^»I^2’ ®'5^' 


T'quation 4.2.8 gives rise 


to the following constraints ; 


S2 S P-i + ^2 ^ ^2 ^ ^ ^2^ ^^4 ^ S ^3 

+ (P, + P2) (P3 ^ S ^4 ^ " ^2) ^^>3 ^ V ^ ^5 

+ SP5 = ® •*’ 

Pj SjS p,-*- p, Sj£ Pj^Si Pp 'h* 

^ h ^2^^5 ^ ^4 ^ + l>4^ ?5 P5 = ° *’ 


= Os 


and 


(p^ , p^) 5 p, . (P^ . P5) S Pa - P, 6 ^5 * " "2^S h 


+ P2 S p ^ “ ^' 


Pol losing 


in^rBolstein's (1) suggestion, let us choose 


5 . , ^*11 .2 
(6 = ^ ^ 
^ i=i 


4 . 2.28 


,. exe ,, i . tl > eposiUo » orthep « leats .- 1 . 
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Therefore, 




a' '^3 


)‘ 


S=1 


4.2.29 


where N(s) is the denominator of z(s) ; 

2 

N(s) = s S 2 + s p^(p^ + p^) + p 2 (p^ + p^) +1 4 , 2.50 

The value of f) for the starting network is 


0 = 2.0 . 

0 

Using this criterion, the optimum shown in fig, 4.2.1(c) is obtained 
The optimum value of 0 is 1.6. It is interesting to note that the 
optimization has resulted in the value E.^ = 0, and thus the 
optimal network contains only as many components as the starting 
canonical realizatiaa. 

The above example brings out an interesting conclusion* 
kf r the criterion involving pole sensitivity, it is possiljs 
that a special member of continu ously equivalent, class, which has 
fewer components than a general member of that class, cam be more 
insensitive than a general member. Leeds and TJgron (IO) have 
conjectured on the basis of experimaital evidence that the 
sensitiv ity decreases as the number of elements increases in 
continuously equivalent network. This conjecture is clearly not 
applicable to continuously equivalent networks in general, as 
demonstrated by the above example. 
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4 •5 Oftiiml Control and Pontryagin' s Principle. 

In this section, we daall discuss tv/o problems of optimal 
control i the terminal control problem and the minimum, integral 
pr'blem. These problems can be transformed into tv/o-point 
boundary value problems by a principle based on variational 
calculus, called the Pontryagin' s Maximum Principle. The Principle 
and its application to the terminal control problem will be 
discussed. A method of converting the minimum-integral problem 
into a terminal-control problem ’/dll be discussed to ^ow the 
generality of the latter problem. A detailed discussion of this 
material can be found in (18). 

Terminal Control Problem 

Given: a system described by the equations 
o 

X = p(x, U) 4*5.1 

^here X is the state-vector of dimension n and U is the input or 
control vector of dimension m; 
the initial state X^ | 
a final time T ; 

and a constant vector B of dimension n, to minimize 

J = BV X(t) ' 4»3*2 

( There pri me denotes transpose) 
by choc sing the apiropriate control 

u® = u®(t) ; t <:t 1.5.3 
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I'ontryagin’ s Maximam Principle : 

The optimal inpat u°(t) minimizing the criterion of 

Iq, 4.5.2 also maximizes the Hamiltonian H(t) given hy 


4.5.4 


at each 


H(t) = P'(t) I'(X(t);lJ(t)) 

• ^ ^ f time 0 t T ; ^‘hexe P(t) is the solution 

ach instant of time u 1 . 


to the differential equations 

iii , 

dt "Sii 

can written tor ccnTeniaace as 

dP il . 
dt ' ’ 

^•ith tie terminal conditions 

4.5.6 

p(t) = 

. jr. rtnYifrol input function is 

Thus, the pxDtlem of searching ora listing 

^ nralue probloiD. consisting 

to solving a two-point Poondary-value p 

„£ Iqs. 4.5.1. 4 . 5 .<. 

ijinimm-integral Problem 

Gi«a the system and the conditions 


4.3.6 


to minimize 


This problem 


'^0 

can he converted into the 


4.5-7 


first problem as follows 


Augmait 


the systaa e 


quations by the equatioa- 



7.ith the initial conditions 

Clearly, 

Tbxis, if -we let 

0 

j • 

b = i * ’ 

1 • 

0 

I (n times) 

1 

L 1 

then for the au;?:mented state equations we have tte first problem. 
Component-level procedures for handling the terminal- control 
problem will be discussed in the following sections. The mininum- 
integral problem has beaa discussed here to an example where 

the first problem covers a frequently used objective, function. 

More such examples can be found in (11, 18). 

4.4 Optimal Inputs in the Time-domain 

Problems of the choice of optimal inputs to dynamical 
systems arise in every field of engineering, from the operation 
of a obemical plant to the directing of a missile . 

'The task of finding the optimal input function, uaially 
from a continuum of permissible inputs, is ' reduced by Pontryagin* s 
maximum principle to that of solving a 2 -point bound aiy-value 
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problem. The maximam principle is especially significant for 
lin .sr systems, since in this case it constitutes a necessary and 
sufficient condition for optimality. 

In the case of linear systems, several advantages of 
retaining the structural details of tte system in the analysis 
have beai shown (6). Among them are the freedom of interconnections 
of subsystems and a physical meaning for each system variable. 

V!e shall outline a procedure for forming, tte state equations 
for an RLC Net orl-:, and show that the corresponding co-state 
equations belong to a similar network, f/hile electrical terminology 
is used for convenience, the analysis is applicable to similar 
non-electrical systems. 

’"e recall from the previous section that the co-state 
eouatiens for a system satisfying the state equations s 
0 

X = F(X,II, t) 44.1 


£ire 



4 . 4-2 


H = P' F ■ 4 •4 -5 

Linear Networks i From the analysis of Chapter 2, we have the 

state eqtaatiDns for linear networks : 



fh' 


rc-^ 

“ 1 ) 

r° 

”■^22 


" 0 

-Ap^ 1 

1 ^5 

! 

( 

[0 

L-' 1 

.4 


0 

+ 

L"®25 

! 

0 1 

j 


0 ir 0 kii- fp olio 


X { U + 5 I j 51 V 

V .[0 V ^®13 °-l 


c"^ 0 T'.r 0 -A. 


i -1 -^ i 

l.o L oj 0 


32 r 


0 G, ! I -B.„ 0 It 

4-1 1 - 12 Jju-^ 5 . 


4 

'n 

J 

0l ■ 
5 

1 B + 

1 

^4!. 


"3 ° ir ° '^ 3 i V, 


0 O^JL-B,, Q 


The co-state eq,uations of a linear system 
0 

X = P,X-QIF 
for linear criteria are 


44.4 


4 . 4*5 


P = - P,'P 

Ti'hich in this case heconie 


0 -A, 


r ®22 


0 -A, 


i-®25 ° 


.R3 ^ 0 0 l 


° -^ 4 ^ Pl 5 " 


r ^ Q |0 - A . 




4 4 *6 


0 
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Substituting the transformation 



44.7 


into Eq . 4 • 5 *6 we get 


d 

dt 

* 

, 

sr 

'c'^ 0 ' 

_ 0 L-'_ 

i)r° 

i '^4 

■ 

+ 

4^ 

1 

0 

1 

1 

'-R_ 0 ' 

i 

i 

1 * 

1 


: [1-^22 V 


0 

K^ 

<M 


L® -«4' 


-1 


i ° ^51 


^ I 


-*13 


■) 

1 

0 

1 

” 0 

/ 

CM 

f 

1 

f" 


J] 1 

° -®4- 


0 J 

1 

. •' I 

1*5- 


44.8 


Equations 4 4*8 resemble Eq, 4 •44* The networks corresponding 
to them are identical eacepfc for two differences t 

(i) All resistance values are made negative 

(ii) ’ Voltage and current driver terms are set to zero. 

7/e shall call the network corresponding to Bq; 44*8 
the adjoint network. It also follows from transfoxmatioi 44 * 7 . 
that the co-state vector P consists of charge and flux variables 
of the adjoint network. On the other hand, if the state variables 
of theoriginal network are chosen to be charge and flux i^ariables, 
then the CO rre spend iig co-state vector consists of the *7oltage and 
current variables of the adjoint network. 

Iifultiterminal Components 

It is best suited to ttiis analysis to treat multi terminal 
com'"onents as they are, that is, to use Eq. 2 .W.3 for their 



terminal equations . In that case the matrix 
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appearing in the above analysis ceases to be diagonal. It can be 
easily shown tint the corresponding matrix in the co-state equations 
is the negative transrose of the original matrix. Thus, for example , 
if a transistor with the terminal equations 



4 . 4-9 


is present in the original net oik, its counterpart in the adjoint 
has the terminal equations 



This result when translated into the language of dependent drivers 
reads as follows. 

Heplace a dependent driver of the oiiginal network with one 
in which the roles of the controlling and controlled edges are 
interchanged, and the parameter is made negative. , 

The type of the dependent driver, charges according to 


table 4 .4 .1 
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Table ■ Typos of dependent drivers in original 

and CO- state netvforks. 


Old Type 

New Type 

■ 1 

4 

2 

2 

5 

3 

4 

1 


fiolution 

In general the optimal control problem requires the 
solution of a t’xo-point boundary value problem, because if for 
the system equations the initial conditions are specified, for the 
ad. joint equations the final conditions are specified. In the 
above case, ho-'ever, the co-state equations are not coupled with 
the state equations, as can be seen from tho fact that the adjoint 
net’ ork is not coupled to the original network .Hence , the co -state 
equations can be integrated backwards in time, from T to 0, starting 
from the value of P(t) given by Eq. 4.5.6 . Then, using the solution 
P(t), the optimal inputs IF° which maximize H can be found and used 
in the state equations yielding the optimal solution of the sy stum. 

It may be noted here that integrating, backwards in time 
is equivalent to integrating forward after having chsnged the sigxis 
of the dynamic components of the adjoint net, voik, That is, the 
signs of components L and C are op-osite to their signs in the 


0 rig Inal net^work . 



other Criteria 
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It ms shovn in Secticn 4*5 that these problems in which 


an integral eriterion is to be minimized can be converted to tm 
gene -all ze--^ mode by suitably augmenting the system equations. The 
co-state equations for the augmented set t 


X = P X + QU 
0 

X = f .(x, u, t) 
n+1 n-H 


4.4.11 


are 


o 

P 


- P'P + s±i 

X 



= 0 


44.12 


Here, the notation 


^ ^n+1 

"Sx 

is used to represent the vector 

riin±i1 

■ < 

! Sx. 

! : 

I 4.. . . 

. * ■ 

^ ^n+1 

in a crncise form. Such a notation in future should be interpreted 


in this manner 
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In ii-q , 4 *4 •'52, P is still an n*-diiiiensionB,l vector* 

The network corresponding to the n-dimensional part of 
4 *5 *12 can he generated from the adjoint network described 


earlier by connecting a current driver of value 


f 

n+1 




c. 

1 


44-15 


in parallel *uth each capacitor and a voltage driver of value 


^ ^n-f1 


4 14 


1 

in series with, each inductor i'hfi orientation of thij drivers 
are sh ."wn in fig, ,4 4-1. 

]i;xainple 44 *1 


In the netrork of fig. 4..4.2(a) find the waveform of the 
driving voltage u(t) which vdll produce the maximum possible 
volt.age across the capacitor C in 10 seconds, given 


i ii('t) ! 4 1 for all t. 

ihe state equations of the above network are 




'0 1 


[''4 


■ o' 


s= 




+ 


i 


-1 -1 .2 

i 



1 

L L- 



. L j 


-1 


’•^'he adjoint network is shorn in fig. 


equations 

in terms of q*^ 

and sue 


^«*c‘ 


' 0 1" 


1*0 ‘ 


>L.^ 


__-1 1.2- 


- K- 


u(t) 4.-V.15 

4.5.2(b) and its 


4.1.16 
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(a) 

The Original network. 


i* = 0.0 


- 1.2 1.0 



(b) 


'^'he adjoint network. 


Pig. 4.4,2 A simple Optimal design problem. 
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'fe arc to minimize - Vp(lO) 


, or 


L-i 0 


( 10 ) 


4.4.17 


therefore 


-1 

0 1 


4.4 .18 


and hence 






( 10 ) = 


1 

0 


4.1.19 


Solving lq.s. d.4.16 tjid 4. 4. 19 we get 


‘n' 

.6(t - 10) 

= 1 .25 6 

mm"] — 

sin (tan 4/2 +8 - .^t) 

[I^L, 


sin (8 - .8t) 

■ _ 


4.1.20 

O / V 

u (,tj is chosen to maximize 


H = 0*^ u(t) 
from which vte get 

u°(t) = sgn (sin (8 - .8t)) 4. 4. 21 

IT sing this value, v^ (IO) (max) is obtained to be 1 .2 volts, 

1 the functi^m signum (sgn) is defined by 

sgn (x) = -1 if X 0 
^ if X = 0 

and = 1 if x 0 



Nonlinear Net wiics ; 
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■ie restate the nonlinear component terminal equations 
from Chapter 2 here for convenience. The mixed form of the 
resistor equations is used for convenience. \'ie have ; 



4.4.22 


for the dynamic components 5 and 



4 . 4 . 23 ' 


for the resistors. If multiterminal resistors are present, their 
terminal equations are included in Eq. 4*4*25* First, we shall 
consider networks of two-terminal components. From Fontryagin's 
principle v!e have the co-state equations : 
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, (say), 44*25 


M = - 


■ d ^*2 




3^2 


V I' 

L 2 5 


B 22 0 


+ I Y' I' I 
L 5 4-* 


"0 

A, ~ 

; 

52 


.-"12 

0 

/ 


4 .4.26 


2144141 .0 

slj 

Bif feraatiating Eqs. 4 *4 •22, 4* 'U23 partially with respect to 


4.4.27 


I ^5 J 


and substituting into Eq. 4 *4 *26 we get 
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0 

“A^ ” 
22 


'■ 0 -A 1 
21 


0 


- "®25 ° - 


U + (-D® ) 
\ r 


-1 


A, 


13 


31 

0 



4 .4.28 


v/hare the meaning of D® and D® has been explained in chapter 2 . 

<e observe that Eq.. 4 . 4.26 is similar to Eq,. 4 . 4,8 where 
the small-signal matrices D® and D® replace the linear component-value 
matrices 


1 -1 



0 

0 


R, Ol 


and 

3 

0 l"'' 


0 G, 

- 


4j 


Therefore, the steps for transforming tte original network into the 
co-state network remain the same. If multiterminal components are 
included, the discussion in the linear case for their replacements 


in the co-state network holds true here also, when small-signal 
characteristics are considered. Therefore, we can give the follomng 
algorithm for converting the original network into the co-state network 
Algorithm 4 .1 .1 . 

(i) Replace all independent voltage drivers by short-circuits and all 
indepc?ndfcnt current drivers by open-circuits. 

(ii) Replace all two-terminal resistors by linear resistors whose 
values are negative of the small-signal values of the original resistors 

(iii) Reverse the roles of controlling and controlled branches for 
all dependent drivers and replace them by linear dependent drivers 
whose value .is the negative of tie small-signal value of original 
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driver and whose type is dictated by Table 4.4. 1 . 

(iv) Replace each capacitor by a linear capacitor vdiose value is 
ecjuol to the small-signal value of the original capacitor, and 
connect in parallel with it a current driver of value 



C. 

1 


'"n +1 




^C. 


1 

wbore is the small-signal value of the capacitor. 

Similarly, replace each inductor by a linear inductor 
whose value is equal to the small-signal value of the original 
inductor, and connect in series with it a voltage driver of value 


L. 

1 


^n+1 


dA. 

1 


where is the small-signal value of the inductor t 
4.5 Tests of Optimality ; 

In the last section the network corresponding to co-state 
aquations was derived from the original network. In the case of 
nonlinear net'-orks the adjoint network depended on ths response of 
the original network, and hence the solution of the two was essentially 
coupled. In the linear case, if the criterion for optimality was also 


linear in tbs fin.3l state, then the two netiTOrks are uncoupled and 
the optima] control problem can be solved as an initial value problem. 

If the problem is known to have a solution, ard. the 
equations have a unique solution, the probltm is completely solved, 
since the obtained solution can only be the optimal one. If eilher 



of these condit ions is- wi£tilfilLed,howerei, there still remains the 
prohlem of determining whether a given solution is optimal or not. 
lot us re-urite the state and co-state equations of the linear case, 
for completeness ; 


o 


X = P X + Q R 

4.5. l(a) 

0 



4.5.1(b) 

? .-P.P-P , 

4.5.2 


■ffe shall discuss some sufficiency conditions for 
optimality for the linear case in order to investigate their possibl 
relcti'n to the topology of the network. The tests discussed belo’w 
are due to Robbins ( 21 ). 

Sufficient conditions for optimality are available for 
the case where- the control vt-jctcr U° to be tested is in the interior 
of the permissible regicn Ui These conditions can, however, be 
applied to a wider class of situations by 'freezing' some of the 
control variabl 8 S or their combinations (See ( 2 l)). 

In the following analysis, we assume that the c aatrol 
vector is in the interior of its region. 

The first condition is that the matrix , 


L ^ u. u. J 


4 . 5-5 


should be negative-definite at 
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O 

for u to be optimal* 

In the system descrihed hy 4.5 -1 this -Corresponds 
to the condition that the matrix 




n-+1 


L 




4-5.4 


he positive definite. If the matrix 4.5.4 is positive non-defini te , 
th,o solution is non-optimal. If it is positive semi-definite, further 
tests must be carried out. 

A further test consists in reformulating the problem 
r.ith possibly fewer control variables so that, for the ne-w problem, 


H 


uu 


= 0 


4.5.5 


(see (I8)), 2'his condition may also arise directly in tie original 

pr-'bl em if the criterion f .docs not de-pend explicitly onlJ, that 

n+1 

is 


f , = f (X, t) 
n+1 n+1 


4 -5 .6 


In the case ’-here Fo . 4.4.5 holds, the following two conditions are 
S’afficient for optimality s 


( AS)) = 0 


4.5.7 


and 


dt^ * 


4.5.8 


is pos it ive def inite . 



ICO 


The condition of Iq, 4 .5 .7 is automatically satisfied by the system 
desor bed by Eq. i.5.1, and the condition of ^q. 4.5 .8 reduces 
to the condition that the matrix 


r 


f 


n+1 




4.5.9 


be positive definite. 

If the matrix is positive non-definite the solution is 


ncn-optiml . If th2 matrix is positive semi-definite, this test 

is also indecisive, for instance when not all the state variables 

influence f . t 
n+1 


t) 5 k n 4.5.10 

B\irther tests can be ap'lied to tMs case; but ife shall not go into 
those details, because they are not necessary for the developaent 
of this thesis. 

In the above discussion we observe that both results 
depend only on the cost hinction, and not ex^’llcitly on the topology 
of the netwf'rk. Thus, it is useless to attempt to give a network 
interpretation for these tests. 

4.6 Further Remarks : 

yie have discussed soma aspects of design which indicate 
compit er ap - licati cn . In particular, parameter optimization can be 
implemented on a digital computer . Some discussion of this topic 
can be found in (12). 
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The advantages of a. network interpretation of co-state 
equations are ma,ny. Apart from all usefulness of the formulation 
as an analysis tool, the procedure has a strong pedagogi.e advantage i 
it mr.kes it possible to explain important results of optimal 
control theory through networic analysis. 

The simplicity of formulation, and its systematic 
nature, mdees both analytical and mechanical formulations of co- 
state equatinns easy and less prove to errors. 

The twj-point boundary value problem in which this 
formulation results can in some restricted cases be solved 
directly by network analysis programnes, as pointed out in section 
4« In some more general cases, it may be possible to solve it . 
iteratively as a succession of initial value problems. 



5. Conclusi'ns : 
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The technology of network construction has developed 
in all facets. The rapid reduction in component sizes and costs 
have made extremely sophisticated electronic equipment aich as 
rocl.et guidance systems feasible. Integrated circuits can be 
ccnsidered as the pinnacle of progress in this direction. Their ■ 
advv^nt caused an immediate change in the philosophy and the means 
of netwoik designers. The need for economy in the number of componints 
is replaced by the need for reliability end the need for design 
•aid analysis of extremely large networks. The requirement of speed 
in analysis has been so great that time-sharing consoles with 
expensive light-pen input and graphical display units are maie 
available at several places (See (3))* 

As design methods better and better suited to digital 
computers are developed, the cost of analysis in terms of computer 
equipm.mt and time is likely to decrease. The present-day design 
oriented analysis progranmes lay emphasis on user interaction. If 
some features of design are made automatic, the need for such 
expensive int'craction v.ill decline. As mentioned in the 
introfucticn the models that have been presented in chapters 3 
and 4 provided the theoretical basis for developing programmes 
\)iiicb would enable to minimize human interaction with the compsitar 
in the design process. 

The foriailation of the methods in those chapters has 
always been kept at the component level. Some advantages of the 
samt; in reference to parameter optimization were pointed out in the 
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conclusi'-n of chapter 4* These advantages, however, have a more 
or less gen'?ral nature. . 

Firstly, computati'rn is systematized by this procedure, sines once 
the transformed network is obtained analysis takes tho same course 
as usual network analysis. ■ . 

The other advantages are all derived from the fact 
that component -level formulation makes possible a physical 
interpretation of the quantities on which com_^utations are made,Thus, 
interpretation of results and error analysis become simpler, aid 
an insight into the role of each component in the problem oai be 
obtained simply by studying the results. 

It can be seen that while the procedures for sensitivity 
and co-state formulations are quite eliigant for linear netvroiks, they 
seam artificial for nonlinear networks, in that they require 
additional information about the network diich is not needed for 
simple analysis. The situation here resembles the attempts to 
modify transform methods in order to make them applicable to non- 
linecr systems. These methods had gained considerable momenhim 
before time-domain analysis became well-known. Likewise, the 
stiudies '^n sensitivi ty and optimization for nonlinear systems will 
be more natural ’'hen new methods are developed than in the present 


studies 
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